Accelerated expansion from a non-minimal gravitational coupling to matter 
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It is shown that a non-minimal coupling between the scalar curvature and the matter Lagrangian 
density may account for the accelerated expansion of the Universe and provide, through mimicking, 
for a viable unification of dark energy and dark matter. An analytical exploration is first performed, 
and a numerical study is then used to validate the obtained results. The encountered scenario allows 
for a better grasp of the proposed mechanism, and sets up the discussion for improvements that can 
lead to a complete agreement with the observational data. 
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I. INTRODUCTION 

Many modifications of theories of gravity are moti- 
vated by one of the outstanding puzzles of modern cos- 
mology: the origin of the observed accelerated expansion 
of the Universe. The most common approach to this issue 
relies on the presence of a dominating dark energy com- 
ponent (with £Ide ~ 70%) [1], which might arise from 
several competing candidates: a cosmological constant 
term in the Einstcin-Hilbert action, a scalar field, usu- 
ally referred to as quintessence @, chameleon fields Q, 
or other alternatives to General Relativity (GR), such as 
the cardassian model with a modification of the Fried- 
mann equation [|| , braneworld scenarios @ or the gener- 
alized Chaplygin gas unification of dark energy and dark 
matter [6|. 

Aiming at a description of this dark energy compo- 
nent, several authors have put forward proposals based 
upon the so-called f(R) models, where a modified ac- 
tion functional exhibiting a non-linear function of the 
scalar curvature R is considered This is usually re- 
garded as stemming from a low-energy phenomenological 
approximation to some higher energy fundamental the- 
ory; indeed, one-loop renormalization of GR requires the 
introduction of higher order terms in the curvature in 
the Einstein-Hilbcrt action functional, and other avail- 
able invariants — such as contractions of the Ricci or 
of the Ricmann tensor — may also arise when quantum 
corrections arising from string theory are considered (see 
Ref. fg| for a thorough discussion). In a cosmological 
context, these models usually rely on a decreasing f(R) 
function that, since the scalar curvature is decreasing, de- 
viates strongly from GR at late times — thus producing 
the required accelerated expansion |{|. 
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Aiming to further extend the f(R) theories, a model 
was advanced exhibiting not only a non-linear f(R) term 
in the action functional, but also a non-minimal cou- 
pling between the matter Lagrangian density C m and 
the scalar curvature [l(|. The purpose of this work is to 
show that this latter model may be used to account for 
the accelerated expansion of the Universe without any 
explicit additional matter component (e.g. scalar fields). 
By resorting to a previous study where it was shown that 
this non-minimal gravitational coupling with matter can 
mimic known dark matter profiles (thus producing the 
reported flattening galaxy rotation curves) , one con- 
cludes that the proposed model yields, through its grav- 
itational impact, a viable unified scheme to mimic the 
presence of both dark energy and dark matter. 

This work is organized as follows: the non-minimal 
gravitational coupling model is discussed in Section II; an 
analytical work, establishing quantitative results related 
to the accelerated expansion of the Universe, is discussed 
in Section III; this discussion sets up the numerical cal- 
culation that confirms the proposed scenario in Section 
IV. Finally, the conclusions are presented in Section V. 



II. THE MODEL 

Following the discussion of the previous section, one 
postulates the following action for the theory [Tpj : 



S = 



- 2 hiR) 



[l + A/ 2 (i?)]£, 



(1) 
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where fi(R) (with i = 1,2) are arbitrary functions of the 
scalar curvature R, C m is the Lagrangian density of mat- 
ter and g is the metric determinant. The contribution of 
the non-minimal coupling of ji is gauged through the 
coupling constant A, which has dimensions [A] = \f2\~ 1 - 
The standard Einstein-Hilbert action is recovered by tak- 



ing f 2 = and /i = 2k(R-2A), 
A is the cosmological constant. 



where k = c jVoixG and 



2 



Variation with respect to the metric g^ v yields the field 
equations, here arranged as 



(Ft + 2\F 2 C m ) Rfj, u — -yfxQixv — 
A^ (Fx + 2\F 2 C m ) + (1 + A/ 2 ) T M „ , 



(2) 



where one defines A^ v = V M Vr/ — g^\3 for convenience, 
and writes Fi(R) = fl(R), omitting the argument. The 
matter energy-momentum tensor is, as usually, defined 
by 



2 S(^C m ) 



V=9 W 
By taking the trace of Eq. ©, one obtains 



(3) 



{F 1 + 2XF 2 C rn )R-2f 1 = (4) 
-3d(Fi + 2XF 2 C m ) + (1 + A/ 2 ) T . 

The Bianchi identities, V^G^ — imply the non- 
(covariant) conservation law 



Guv 



1 



(8) 



so that, using L r , 
defines 



2k Miy p ' 
-p (see [l3[ for a discussion), one 



T, 



2k 



T c 
(11/ 



Ft - 2F 2 p 
2k 

Ft - 2F 2 p 
1 



(9) 



A fil/ (F 1 -2F 2 p) 



-(ft - FtR)g^ + F 2P R 9llv + / 2 T„„ 



One assumes that the matter content of the Universe 
is described by a perfect fluid, endowed with an energy- 
momentum tensor 



Tfu, = (p + + pg^ 



(10) 



Resorting to the Friedmann-Robertson- Walker metric 
given by the line element below, 



XFo 

V% = t- ±r(gp,£ m -Tp,)Vi*R , (5) 

1 + AJ2 

which, in the context of an analogy between Eq. ([T]) 
and a scalar-tensor theory, may be interpreted as due to 
an energy exchange between matter and the scalar fields 
associated with the model [l[ for the non-trivial ft(R) 
and f 2 (R) terms [13 • 

Since a complete study of the joint effect of a non- 
trivial ft(R) and f 2 (R) is too involved, one focus the 
attention on the latter, thus setting ft(R) = 2kR (dis- 
carding the cosmological constant A); this reduces Eq. 
© to 



1 + -F 2 C m ) R^ - \Rg, lv = (6) 



- (CW - 9^U) (F 2 C m ) + -!- (1 + A/ 2 ) T M „ 
and, taking the trace, the equivalent of Eq. (gj), 



1 - -F 2 C m )R = (7) 

K 



3-D(F 2 £ m )--^(l + A/ 2 )T 



III. ACCELERATED EXPANSION PHASE 
A. Power-law expansion 

One begins by rewriting Eq. ^ in a more natural way, 



ds 2 = -dt 2 



a 2 (t) 



dr 2 



VI - kr 2 



(11) 



and considering the isotropy and homogeneity implied 
by the Cosmological Principle, one writes the four- 
velocity as — (uo, 0,0,0). The normalization con- 
dition u^u^ = — 1 thus translates into Uq = 1, and the 
energy-momentum components reads 



Too = P j 



(12) 



a p 



It is easy to check that, likewise T M „ with the adopted 
four-velocity, T° is also diagonal: one may write the 
corresponding curvature "density" and "pressure" , 



2k 



(13) 



(h-F 2 R)p--(h-FtR) 



Ft - 2F 2 p 
3H \(F[ - 2F 2 p)R - 2F 2 p 



T 



Pc 



(14) 



2 k 



(F{-2F 2 , p)(R + 2HR) 



Ft - 2F 2 p 
(F{' - 2F' 2 'p)R 2 + F 2 [pR - 2(p + 2Hp)] 

l -(ft-F 1 R)-^F! 2 pR + f 2 p 
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Defining the Hubble parameter H = a /a, leads to the 
Friedmann equation 



H 



k 1 , 



6k 

and the Raychaudhuri equation 



(15) 



a 



1 



= H + H< = -—[p n 
a 12k 



+ Pc + 3(Pm + Pc)} , (16) 



the latter providing the condition for an accelerated ex- 
pansion of the Universe, p m + p c + 3(p m + p c ) < . 
Defining the deceleration parameter as 



aa 

q = -T7T , 



(17) 



one may combine the two preceding equations to obtain 



]_Pc_ 

4k H 2 



(18) 



To provide an insight on the effect of the coupling with 
matter, one may write the curvature pressure and density 
p c and p c arising from two separate cases. On one hand, 
in the absence of matter, f%{R) — 0, one has 



Pc 
Pc 



2k 
F~i 



-(f 1 -F 1 R) + 3HF{R 



(19) 



2k 



F[(R + 2HR) + F['R 2 + - Frf) 



On the other hand, including matter, f 2 (R) ^ 0, and 
setting fi(R) = 2kR, leads to 



K 



Pc 
Pc 



k - F 2 p 

K 



(f 2 -F 2 R)p + 6H [F! 2 pR + F 2 p 
- 2F 2 p(R + 2HR) - 2F 2 pR 2 + 



k — F 2 p 

F 2 [pR - 2{p + 2Hp)} - pi? + f 2 p 



There is a clear difference between the expressions above: 
the latter two explicitly depend upon the matter density 
and pressure p and p, while the former two are a func- 
tion of R (and its derivatives) only. This, of course, stems 
from the non-minimal coupling between matter and ge- 
ometry, but implies that one cannot simply neglect the 
contribution from matter when solving the Friedmann 
and Raychaudhuri equations. Indeed, in Ref. [f|, this 
allows for the determination of a relation between the 
evolution of the scale factor a(t) and the exponent m 
present in the non-trivial term fi(R) — R\(R/ Ri) m \ 
in the present case, although one can still assume that 



Pm < Pc P < \Pc\, the density and pressure appear in the 
above definitions, so that this does not translate directly 
into setting p = p = 0. 

In order to solve Eq. (| 15[) and (|16p , one assumes a flat 
k = scenario and inserts the Ansatz a{t) = <Xo(t/to) for 
the evolution of the scale factor — physically interesting 
since it gives rise to a constant deceleration parameter, 
with (3 > for an expanding Universe and /3 > 1 for 
accelerated expansion. Thus, 



H 
R 



a 
a 



t 



(21) 



aa 

• 2 



i 



a 
a 

- 1 



6(H + 2H 2 



Since one is interested in studying the effect of the 
non- minimal coupling f 2 (R), one assumes that it domi- 
nates any non-trivial addition to the usual linear curva- 
ture term, thus one can set for simplicity fi(R) = 2kR. 
Given the assumed power-law expansion, the former is 
assumed to have the form f 2 (R) = (R/R 2 ) n , prompting 
for the search of the relation between the exponents n 
and /?, as well as the physical meaning of the coupling 
strength R 2 . 

Furthermore, one requires some foreknowledge of the 
evolution of the matter density p and pressure p; one 
may assume that it is modelled as a dust distribution 
with p — 0, and looks at Eq. ([5]) for the evolution of 
p. Luckily, although the non-(covariant) conservation of 
the energy-momentum tensor is perhaps the most strik- 
ing fundamental implication of the model here studied, 
the introduction of the adopted energy-momentum for a 
perfect fluid and the Lagrangian density C rn — — p yields 
a vanishing r.h.s. for Eq. ([5]). 

Following the standard interpretation, one can state 
that the expansion of the Universe remains adiabatic, 
with no direct transfer of energy between matter and 
"curvature" component expressed in Eq. ([20| . The v = 
component of Eq. (0 reads: 



p + 3Hp = 0^ p(t) = p 



o(t) 



t 

P0[J 



3rf 



(22) 



With the above law for p(t) and the expressions for 
f 2 (R) and a(t), the curvature density p c and pressure p c 
read 



p c = 6kpq/3 x 



(23) 



1 - 2(5 + n(50 + 2n - 3) 



(&) *-(&) (£) [W-i)i IH, « 



and 
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p c = 2Kp n x 



(24) 



2 + An 2 - 0(2 + 3/3) + n(80 - 6) 
^-(i) (i) [6/3(2/3 -l)] 1 ""* 

— 1/2 

defining t2 = R 2 , for simplicity. 

In what follows, one assumes alternatively that the de- 
nominator of the curvature pressure and density is domi- 
nated by either F 2 p or the constant k; for simplicity, the 
regime F 2 p > k is dubbed "+ regime", with the converse 
leading to the "— regime" . Actually, one sees that 



-2p n 



2 + 4n 2 - 0(2 + 3/3) + n(80 - 6) 



Since the l.h.s. of the Fricdmann Eq. (|T5|) falls as i~ 2 , the 
exponent can be directly obtained from the expression 
for the curvature density, 



3/3 + 2n : 



= 0_(n) = -(l-n) 



(28) 



Inserting this back into the Friedmann equation yields 
the value for the initial density, 



F 2 p 



R 
R~ 2 



P_ 

R 



n [6/3(2/3 -l)]"- 1 ^ 



t zn f 







£2(n-l)+3£ 



thus scaling as t 



2-2n-3/3 



(25) 



The + regime: Fip > ft 



Clearly, both the curvature pressure p c as well as the 
density p c experience two separate time evolutions, sig- 
naled by the relevance of F 2 (R)p on the denominator of 
Eq. (HUD. 

One first attempts to solve the Friedmann and Ray- 
chaudhuri Eqs. (fT5|) and (fT6|) in the regime F 2 p > k, 
so that the curvature pressure and density are given ap- 
proximately by 



6k0 (1-20 
Pc = ~jT ( + 5/3 + 2n - 3 



(26) 



2 k 



Pc 



— [2 + 4n 2 - 0(2 + 3/3) + n( 



6) 



Inserting this into the Friedmann Eq. (|15l) leads to = 
+ = (1— n)/2, which trivially satisfies the Raychaudhuri 
Eq. (|16p . The condition for an expanding Universe 0+ > 
yields the constraint n < 1. 

Replacing onto Eq. (l25t leads to F 2 p oc since 
0+ > 0, one concludes that, once the inequality F 2 p > k 
sets in, the l.h.s. increases with time: the + regime, once 
attained, remains valid. 



2. The — regime: F^p < a 

The analysis of regime F 2 p < k is slightly lengthier, 
since the value of po appears explicitly in the approxi- 
mated expressions for curvature density and pressure, as 
can be seen below: 



-Gpo0- 



l-20 + n(50 + 2n-3) 



(27) 



/- 



(29) 



As before, the above expressions trivially satisfy the Ray- 
chaudhuri Eq. flTB"]) . The condition for an expanding 
Universe 0- > also leads to the upper bound n < 1. 

Since 3/3_ = 2(1 - n), Eq. (gS]) indicates that F 2 p is 
constant. Therefore, one concludes that if the — regime 
is attained it is also permanent. 



3. Regime validity 

In the previous paragraphs, one has concluded that 
there are two possible regimes where the effect of the 
non-minimal coupling f 2 is dominant, corresponding to 
the positive or negative sign of F 2 p — k. Furthermore, 
it was shown that the + regime leads to an increasing 
F 2 p term, so that the corresponding inequality F 2 p > k 
becomes even stronger; likewise, the — regime yields a 
constant F 2 p term, so that F 2 p remains smaller than k. 

If there are no possible transitions between the two 
regimes, there is still the issue of the onset of the dom- 
inance of the non-minimal coupling, that is, the condi- 
tion f2{R) > 1- Since one adopts the power- law form 
HiR) = {R/R2) n and the scalar curvature R decreases 
and one aims for a late-time dominance leading to the 
currently observed accelerated expansion, one concludes 
that the exponent n must be negative. 

This requirement for a negative exponent is, in essence, 
analog to the one found in a previous study concerning 
a mimicking mechanism for dark matter [ll| : in that 
work, an inverse power-law was required so that the ef- 
fects of the non-minimal coupling become dominant at 
large distances, thus leading to the flattening of the rota- 
tion curves of galaxies. Conversely, the application of the 
considered model to astrophysical objects with high den- 
sities (such as the Sun [l4j ) leads one to the consideration 
of a linear coupling f2(R) oc R. 

Since R = 60(20 — l)/i 2 , it is clear that the non- 
minimal coupling dominates after a transition time T = 
y/60(20 — l)t 2 (when R < R 2 ); considering that is of 
order unity (0 = 2/3 for the matter dominated phase), 
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this translates into T ~ t 2 . Thus, the issue of evalu- 
ating the sign of F 2 p — k simplifies to the direct eval- 
uation of this quantity at T = t^- To do this, one re- 
quires the parameters to and po specifying the evolu- 
tion of the matter density: it suffices to consider the 
WMAP7 value to = 13.73 Gy and po — Q m p C rit, with 
~ 0.3 the relative matter density in the Universe 
and p cr it = 3Hq/8itG ~ 10~ 26 kg/m 3 its critical den- 
sity, from Ref. [Hj]; replacing k = l/(16irG) yields 
Potf)/n ~ 0.1, thus 



F 2 p/K 



\F2P\t=U = 



n [6/3(2/3 - I)]"" 1 



pq4 



O.lna 11 - 1 (4n 2 - 5n + l) 



n-1 



3(3-2 



* (30) 

3/3-2 




FIG. 1: Value of F 2 p/n (Eq. 
(black) or F 2 p < k, (gray). 



32])), assuming that F 2 p > K 



where a signals the alternative regimes, a = 1 for the + 
regime or a = 4/3 for the — regime. 

Since /3 > 2 (marking an accelerated expansion of 
the Universe), one concludes that an earlier onset of 
the non-minimal coupling dominance (that is, a smaller 
value of t 2 ) increases the r.h.s. of the above expression 
- eventually leading to the condition F 2 p/ 'k > 1 that 
marks the -I- regime. Indeed, there is an interplay be- 
tween the two terms of the above expression: the term 
0.1n[4(4n 2 — 5n + 1)/3]" _1 is smaller (in absolute value) 
than approximately 5.6 x 10 -3 in the domain n < 1; how- 
ever, since the accelerated expansion has already begun, 
t 2 < t Q and /3 > 1, thus leading to {t /t 2 ) 3l3 ~ 2 > 1. 

The issue is settled by resorting to observational data 
in order to fix the accelerated expansion onset time 
tfi ~ ~ti- This may be written in terms of the redshift ze 
marking the change of sign of the deceleration parameter 
q, through 



3-2/3 



3/3-2 



3-2//3 



= (l + z E ) 



3-2//3 



(31) 

Taking the best-fit value ze ~ 0.36 [l6[ and inserting 
into Eq. (|3T)1) . together with the expressions for a and 
/3_(n), p+(n), yields 



B. Discussion 

The previous section has allowed one to conclude that 
the effect of the non-minimal coupling leads to an evo- 
lution obeying condition F 2 p < K, so that the scale fac- 
tor scales with a power-law a(t) oc t^(n), with /3(n) = 
/3_(n) = 2(1 — n)/3. The assumption of accelerated ex- 
pansion q < —¥ /3 > 1 (more stringent than the previ- 
ously assumed /3 > 1, where no acceleration is required) 
leads to the stronger constraint n < —1/2. One may re- 
sort to the expression for the deceleration parameter in 
Eq. ([22]) and write instead 



(33) 



2(1 + 9) 



2(1 -n) 



Clearly, as n — > -co, one gets q — > — 1: this is the ex- 
pected value if one assumes that the effect of the non- 
minimal coupling mimics the ACDM scenario, thus repli- 
cating a cosmological constant (although that implies an 
.exponential evolution of the scale factor, not a power-law 
one). 

One may easily evaluate the equation of state (EOS) 
parameter w, as given by relation p c — wp c ; manipulat- 
ing the Friedmann Eq. (fT5|) together with Eq. (fT5f yields 
the well-known relation 



-F 2 p = 0.1nUn 2 -5n + l) n 1 x (32) 

K 

r ( i. 36 )(i+3«)/(»-i) , F 2 p > k 

| (L36) 3n/(«-i) (4/ 3 )»-l ,F 2 p< K 

The two curves above are depicted in Fig. [TJ for negative 
exponent n: clearly, both are (in absolute value) always 
below unity, so that the condition F 2 p > k is inconsistent, 
while the converse is valid for all n. Hence, one concludes 
that the — regime F 2 p < k is followed by the system once 
the non-minimal coupling becomes dominant. 



Since the exponent n is negative, one concludes that the 
EOS parameter of the mimicked dark energy obeys — 1 < 
w < 0, thus fulfilling the weak, dominant and null energy 
conditions (for am extended discussion, see Ref. (l7|)- 

As stated in Ref. [11], one may assume that the non- 
minimal coupling comprises several contributions, with 
a power-law referring to the dominant term (possibly of 
a Laurent series of a more evolved form) in a particular 
context. Hence, it is not required that the exponents 
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n of cosmological relevance play a part in astrophysical 
contexts, and vice-versa. 

This said, the applicability of the model here con- 
sidered to the puzzle of the flattening of the rotation 
curves of galaxies was studied in a previous work [Tlj j ; 
in particular, it was found that the Navarro-Frenk- 
White and isothermal dark matter profiles can be derived 
from power-law non-minimal couplings with exponents 
tinfw = — 1/3 and njs = — 1, respectively. 

Notice that the constraint n < —1/2 rules out a 
cosmologically relevant non- minimal coupling f2(R) = 
(i?/i?3) -1 / 3 : this is in agreement with the results ob- 
tained in Ref. [ll|, where it was shown that the char- 
acteristic length scale r 3 = 1/ \fR$ is much smaller than 
the relevant Hubble radius Th- 

Conversely, a simple inverse coupling (R) = Ri/Ris 
allowed and yields an asymptotic deceleration parameter 
q = —1/4. From Ref. jl6j . one finds that this value lies 
still within the 2a interval for the present value of q(t); 
however, in Ref. [ll[ it was found that, although such 
inverse coupling could play a cosmological role, it is not 
a dominant one (that is, the characteristic length scale 

n = i/VW<r„). 

Hence, one assumes that the n = — 1 scenario does not 
correspond to the observed cosmological dynamics, and 
an even smaller value of the exponent n is required (thus 
leading to a larger value of a fully consistent model 
would require that the cosmologically relevant coupling 
(characterized by an exponent nc) does not disturb the 
dark matter mimicking scenario already obtained with 
the aforementioned exponents unfw = — 1/3 and nj$ — 
— 1. This shall be the object of a future study, as it is 
clearly outside the scope of the present work. 



IV. NUMERICAL RESULTS 

In this section, a numerical evaluation of the solution 
to Eq. (TTSJ) is performed, varying both the time scale 
t% = l/y/R~2 and the exponent n. This equation is cho- 
sen, instead of the Friedmann Eq. (|15[) or Raychaudhury 
Eq. (TT6")) for three reasons: it best expresses the relation 
between a negative pressure and an accelerated expan- 
sion; it requires only the evaluation of p c , not p c ; it allows 
for a direct comparison with available data. 

In the context of this work, a natural candidate for an 
observable quantity is the evolution of the deceleration 
parameter with the redshift; since the proposed model of- 
fers a clear mechanism for the transition from the matter 
dominated phase (characterized by q — 1/2) to an ac- 
celerated expansion regime, one aims at comparing the 
solution to Eq. (fl~8| with available q(z) evolution curves 
[l6j (see also Refs. [HI), which employ fitting functions 
exhibiting an asymptotic future behaviour. It is useful 
to notice that these studies indicate that the deceleration 
parameter has not yet evolved completely to the perma- 
nent regime q — ¥ const.] the best fit found for the present 
value of the deceleration parameter go ranges from —0.76 
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FIG. 2: Numerical result for F2p/n, for the first case (n — —4, 
ti = to/4, dashed) and second case (n = —10, t% = to/2, full). 



to absolute values larger than unity. 

With this considerations in mind, one may qualita- 
tively predict the impact of varying the model param- 
eters ti and n: increasing the former shifts the transi- 
tion time tx oc t%, i.e., decreases the transition redshift 
Ze, defined by q{zE) = and found to lie in the range 
0.2 < ze < 0.4. From Eq. ([33|) . decreasing the nega- 
tive exponent n lowers the asymptotic q(n) value for the 
deceleration parameter. 

At first glance, it appears that the problem of finding 
more suitable values for n and ti is not too difficult, since 
one may "guide" these quantities to match the reported 
values for go an d Ze- However, the situation is somewhat 
more evolved: asides from obtaining the desired q{z) pro- 
file (roughly specified by these quantities) , one must also 
verify that the Hubble parameter H and the matter den- 
sity (or, equivalently, the scale factor ait) cx p^ 1 ^ 3 ) ac- 
quire their present values. Furthermore, one should guar- 
antee that the time assigned at the present does not de- 
viate widely from t , i.e. H{t now ) ~ H , p{t now ) ~ p , 
with t now ~ to- 

This said, the prospect of obtaining a completely co- 
herent picture with only the postulated power-law non- 
minimal coupling appears unattainable. For this rea- 
son, one should clearly restate the purpose of this work: 
not to provide a thorough matching to the observational 
scenario, but to describe a possible mechanism through 
which a non-minimal gravitational coupling might ac- 
count for the key features discussed: a transition from 
a matter dominated to a Universe with an asymptotic 
accelerated expansion. For this reason, two numerical 
studies were undertaken, as detailed below. 

The validity of the + regime was also ascertained, with 
F%pl 'k found to lie below the 10 -33 level. As can be seen 
in Fig. [5J this quantity is decreasing, instead of being 
constant, which does not contradict the finding of the 
previous section, since condition F^p = const, was de- 
rived in the context of a constant deceleration parameter, 
a regime which has not been attained yet. 

For definitiveness, in Fig. [5]one compares the solutions 
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q(z) to Eq. §TE\) with the profile given by the fitting of H ^s'lMpc) 
the function 



/ \ 1 Q\Z + q 2 
?(Z)= 2 + (IT? 



(35) 



to the Sloan Digital Sky Survey and the WMAP7 [15| 
combined data, with best fit parameters q\ = 1.47 and 
q 2 = —1.46 [16[ (Fig. 1 therein). The evolution of the 
EOS parameter w, derived from Eq. (|34p . is presented in 
Fig. [6l Both figures show the solutions flowing to nega- 
tive redshifts, thus showing the onset of the asymptotic 
regime. 



First case 




FIG. 3: Evolution of the Hubble parameter H(t) for the first 
case (n = —4, t2 = to/4, full) and the second case (n = —10, 
t2 ~ to/2, dashed). 



This first numerical exploration aims at obtaining the 
evolution profile of the deceleration parameter q(z), con- 
strained by the requirement of matching the matter den- 
sity and the Hubble parameter, to their currently ob- 
served values at a time t now = to (as can be seen in Figs. 

By varying the exponent n and the characteristic 
timescale t 2 , it was found that this is obtained for n = —4 
and t 2 = to/ '4. From Eqs. (f3"3"|) and (IMl) . one sees that 
the exponent n — —4 yields an asymptotic deceleration 
parameter q(n = —4) = —0.7 and an EOS parameter 
w(n = -4) = -0.8. 

From Fig. [5j one sees that the agreement with Ho 
and po is attained at the expense of a higher value for 
<7o = —0.53 and the transition redshift ze = 1- The 
solution q(z) misses the indicated marks qo < —0.76 and 
0.2 < ze < 0.4 and falls mostly within the 3a allowed 
region (except in the vicinity of z — 0.3); the region z < 2 
falls within the la region. 

B. Second case 

The second attempt relaxes the above constraint, and 
aims instead to obtain a transition redshift within the 
range 0.2 < ze < 0.4, with most of the solution falling 
within the la allowed region. This is obtained for 
n = -10 and t 2 = t /2. From Eqs. (J33J) and d34j), this 
value for the exponent yields an asymptotic deceleration 
parameter q(n — —10) — —0.86 and an EOS parameter 
w{n = -10) = -0.91. 

Following the previous discussion, Figs. [5] shows the 
converse trade-off: the transition redshift ze — 0.36 is 
obtained, and the solution q(z) closely approaches the 
best fit curve before it, z > ze\ the current value for qo 
falls within the 2a region. However, both the Hubble pa- 
rameter as well as the matter density deviate from their 
present values, H(t ) = 0.78H and p{to) — l-6po- 

Alternatively, one may express this mismatch by deter- 
mining the times at which these quantities attain their 
present values: it is found that H(0.7to) — Ho and 



Pel Pcrit . Pel Peril 




t/t 



FIG. 4: Evolution of the curvature density p c (black) and 
pressure p c (gray) for the first case (n = —4, t2 = £o/4, full) 
and the second case (n = —10, t2 = to/2, dashed). 

p(1.2to) = po- Clearly, this result does not amount to 
a simple shift of t now with respect to to, and expresses 
the aforementioned over determination of the observable 
quantities with respect to the available model parameters 
n and t 2 . 



V. DE SITTER SOLUTION 

In the previous sections it was assumed that the Uni- 
verse evolves towards an accelerated expansion phase 
due to a power-law non-minimal gravitational coupling 
f 2 (R) = (R/R 2 ) n , tailored from the assumed Ansatz 
for the scale factor a(t) oc t 3/9 . The latter does not 
allow for a De Sitter phase characterized by q = — 1, 
which corresponds to an exponentially evolving scale fac- 
tor a(t) = aoexp(Hot). From Eq. (|33|) it is found that 
this regime is approached in the limit n — > — oo. 

For completeness, one may instead investigate the par- 
ticular form f 2 (R) that gives rise to a de Sitter spacetime. 
Since this yields a constant Hubble parameter H = Ho 
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may proceed and consider that the above equation hints 
that the non- minimal gravitational coupling f2{R) must 
follow the aforementioned "strong" + regime, F 2 p > k. 

It is trivial to argue that this regime cannot occur for- 
ever: the De Sitter phase is characterized by a constant 
scalar curvature vis-a-vis a constant i^(-R); however, the 
matter density p decreases monotonically, so that F 2 p 
also drops. Hence, even if the condition F 2 p > k is ver- 
ified at the onset of the exponential expansion phase, it 
will eventually be untenable. 

One must also check either the Friedmann Eq. (|15[) or 
the Raychaudhuri equation (|16p. for consistency. From 
Eq. ([55]). one has 



FIG. 5: Evolution of the deceleration parameter q(z) for the 
first case (n = —4, ti — to/4, full) and the second case (n = 
—10, ti = to/2, dashed); from Ref. la, 2a and 3a allowed 
regions are shaded, white line gives best fit. 
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FIG. 6: Evolution of the EOS parameter w(z) for the first 
case (n = —4, ti = to/4, full) and the second case [n = —10, 
t2 = to/2, dashed). 



and scalar curvature R = 12Hq, the curvature density 
and pressure (Eqs. (fT5|) and Q3|), respectively) take the 
form 



Pc = 
Pc = 



Kp 



k - F 2 p 
up 

K - F 2 p 



(h - 30tf 2 F 2 ) 
6H*F 2 . 



(36) 



To ascertain the form for the non-gravitational cou- 
pling f 2 (R) , one inserts the above relations into Eq. (fl~8| , 
obtaining 



_ 1 J_ Vc_ 

q 2 + ak m 



-1 = 



F 2 p 
k - F 2 p 



(37) 



which has no exact solution, hinting that no dependence 
f 2 (R) will yield the result q = — 1. Nevertheless, one 



Pc + 'iPc 



Kp 



K - F 2 p 



(/a - 12if 2 F 2 ) ~ k ( 12H 2 - A 



(38) 



considering the previous condition F 2 p > k for the 
regime. 

Hence, the Raychaudhuri equation becomes 



1 + h = 24F 2 ff 2 
which admits the solution 

l + f 2 (R) = Kexp 



R 
~R 2 



(39) 



(40) 



defining R 2 = 24Hq and with K an integration con- 
stant. This is clearly a unphysical result, since K must 
be positive (so to match the Einstein- Hilbert action when 
R <C R 2 ), one obtains an increasing function f 2 (R) of the 
scalar curvature. This becomes less and less relevant as 
R decreases during the matter dominated phase, and the 
transition to the De Sitter phase never occurs — on the 
contrary, the non-minimal coupling dominates at an early 
epoch. 

For illustrative purposes, one entertains the following 
possibility: is it feasible to "push the envelope" in order 
to obtain a different form for f 2 (R), by reinterpreting Eq. 
([31))) ? One may assume that the simple De Sitter phase is 
not actually enforced in nature, and the scalar curvature 
will not be constant: one could perhaps replace the factor 
24H$ -> 2R, thus recasting Eq. ([39]) as 



l + .f 2 = 2RF 2 . (41) 

However, this also yields an increasing solution 1 + 
f2(R) — ±y/R/R2, with R 2 a free parameter. Again, 
one must select the positive, increasing solution, with 
the same unphysical result. 

Clearly, no amount of creativity can go against the sim- 
ple interpretation of Eq. ([39| : since 1 + f 2 must be pos- 
itive, so should F 2 . Hence, no form for the non-minimal 
gravitational coupling f 2 (R) yields a De Sitter phase. 
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In order to close this session, one may resort to a simple 
argument that corroborates this impossibility: the sim- 
plest model that originates a De Sitter phase relies on 
the presence of a cosmological constant. Hence, one may 
naively expect that the non-minimal gravitational cou- 
pling term should behave as a constant, (1 + j2)C m ~ 
const. 

Recycling a previous point, one remarks that the con- 
stant scalar curvature yields a constant term /^(-R), so 
that this mimicking of the cosmological constant would 
demand a non-evolving Lagrangian density C m \ Neither 
the adopted choice C m = —p or any of the classically 
equivalent "on-shell" forms (i.e. C m = p [HI) remain 
constant, since the relevant thermodynamic quantities 
decrease as the Universe expands. 

VI. CONCLUSIONS 

In this work, one has applied a model exhibiting a non- 
minimal gravitational coupling with matter to the funda- 
mental issue of the observed accelerated expansion of the 
Universe. In order to do so, one first assumes a constant 
deceleration parameter — 1 < q < 0: the related power- 
law form a(t) oc t 13 of the scale factor leads one to consider 
a power-law non-minimal coupling /2CR) = {R/ R2)~ n > 
so that a negative exponent is required to drive the tran- 
sition away from the early matter dominated phase. 

In the analytical study, one first derived the equivalent 
form for the Friedmann and Raychaudhuri equations; 
since an exact solution does not exist, the two competing 
regimes F%p > ft or F^p < k were discussed, and it was 
found that the later dominates throughout the evolution 
of the Universe. The identification f3 = 2(1 — n)/3 was 
thus obtained, as well as the dependence of the deceler- 
ation parameter q = — 1 + 3/[2(l — n)] and the related 
EOS parameter w — n/(l — n). 

The obtained scenario was numerically tested through 
the variation of the parameters n and R2 . The available 
evolution profiles for the deceleration parameter q{z) and 
the EOS parameter w(z) were compared with the ob- 
tained solutions, and it was found that a thorough fit 
of these results is not fully compatible with the current 
values for the Hubble parameter and the matter density. 

However, the driving force behind the present effort 
is not to yield a complete match with observations, but 
to thoroughly explore the proposed mechanism leading 
to an asymptotic accelerated expansion Universe. For 
this reason, the simplifying assumptions of a linear form 
for the pure curvature term fi(R) — 2kR was made, 



trading the perceived loss of flexibility that a combined 
analysis of non-trivial /1 (i?) and /2 (R) might yield with 
the advantage of deriving analytical results. 

With this in mind, the present work should be regarded 
as a first step towards a more complete description based 
on the encompassing concept of modifying the Einstein- 
Hilbert action, and in a similar way to what was per- 
formed regarding the flattening of the galaxy rotation 
curves, a future improvement would include the effect of a 
combination of power-laws for /2(i?). More ambitiously, 
one could attempt to reverse-engineer the exercise and 
read the form for fi(R) and f2(R) from the observed 
evolution profiles, using the latter as inputs, instead of 
targets. 

Despite the above justification, one could consider that 
the results that were used for comparison (based on Ref. 
[l|| ) employ a fitting function that might be proven in- 
adequate. Indeed, the main qualitatively difference be- 
tween these scenarios and the solutions obtained in this 
work is that the latter present a much smoother transi- 
tion from the matter dominated to the accelerated ex- 
pansion phase: speculatively, one can state that perhaps 
the relatively small number of well selected supernovae 
observations close to the recognized transition redshift 
ze ~ 0.3 might allow for a shallower q(z) transition. 

A final remark is in order: earlier in the text, it was 
suggested that the versatility of the proposed model nat- 
urally accounts for the effects of the presence of the 
two "dark" components of the Universe: dark energy 
and dark matter. Indeed, the possibility that different 
terms present in the non- minimal coupling f2{R) man- 
ifest themselves at distinct scales, astrophysical or cos- 
mological in nature, does not require that the same ex- 
ponent n and characteristic scale R n is shared between 
widely differing phenomena. 

In the present case, it was shown that the n = — 1 or 
n = —3 exponents, which have been shown to account 
for the observed flattening of the galaxy rotation curves, 
do not have cosmological relevance — allowing for other 
values for n to assume such role. Notwithstanding, it is 
interesting to notice that the same EOS form arises both 
for dark matter as well as for dark energy, given by the 
parameter w = n/(l — n). 

This unification is also illustrated by the formulation of 
the proposed model as a multi-scalar-tensor theory [r|: 
in this context, an interesting target for future research 
lies in the possibility of bridging the mechanism here de- 
scribed with others encountered in the literature, namely 
quintessence models Q, chameleon fields [3[ or, the afore- 
mentioned generalized Chaplygin gas unified model @. 
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